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We investigate the impact of the addition of nanoparticles (NP) on the fragility of a model glass- 
forming polymer melt by molecular dynamics simulations. We find significant changes in fragility for 
nanoparticle volume fractions (j> exceeding « 5 %, where fragility changes correlate with the inverse 
variance of the magnitude of the Debye- Waller factor {u^), a measure of local "stiffness" fluctuations. 
We also conflrm the validity of the Buchenau relationship between {v?) and the structural relaxation 
time T for all (j) and polymer-NP interaction types. 



The addition of a small concentration of nanoparticles 
(NP) to glass-forming polymer materials can lead to large 
property changes that are difficult to comprehend by ex- 
tension of the effects of macroscopic filler additives and 
the corresponding theory of composite materials [H [5]. 
For well-dispersed NP, these changes are often rational- 
ized by the large surface-to-volume ratio of these parti- 
cles, which results in a greater interfacial interaction with 
the surrounding polymer matrix. However, this is not a 
unique mechanism for all property changes, and other 
mechanisms, such as chain bridging [31I31[S] and NP self- 
assembly into extended structures [6], are under current 
discussion. 

Of the many affected properties, changes in the glass 
transition temperature Tg have been particularly empha- 
sized, since changes in Tg are correlated with changes 
in diverse transport phenomena. Specifically, both ex- 
perimental and theoretical studies have indicated a ten- 
dency for highly attractive or repulsive (non-attractive) 
polymer-NP interactions to increase or decrease Tg, re- 
spectively. This phenomenon has been rationalized in 
terms of the influence of the NP boundary interactions 
on the dynamics of polymers within an interfacial layer 
near the NP surfaces [3 [HI [HI HD]- In particular, the 
interfacial polymer layer around the NP shows a slow- 
ing down (increased Tg) or acceleration of dynamics (de- 
creased Tg) when the polymer-NP interactions are attrac- 
tive or repulsive, respectively. These observations led to 
the suggestion that the interparticle distance (related to 
the particle concentration for uniformly distributed NP) 
plays a role analogous to film thickness in thin polymer 
films [30. 

Changes of Tg, while informative, provide only a lim- 
ited understanding of how NP affect the properties of 
glass-forming polymer melts. It is also natural to expect 
that the temperature dependence of the dynamic prop- 
erties approaching Tg, referred to as the "fragihty" of 
glass formation [IT|[T2], will be altered. Fragility changes 
have been argued for on theoretical grounds [T2], based 
on the finding that any factor that influences the molec- 
ular packing in the glass state (T < Tg) should also al- 



ter the fragility of glass formation. The NP we study 
should be particularly effective at modifying molecular 
packing, since their size is roughly commensurate with 
the heterogeneity scale of fluids near their Tg (i.e., 2 nm 
to 3 nm) ^13] . 

The present work addresses how polymer-NP interac- 
tions and NP concentration affect the fragility of glass 
formation, and how fragility changes relate to variations 
in the high frequency molecular dynamics, as measured 
by the Debye- Waller factor {v?). We find that the ad- 
dition of nanoparticles can increase or decrease fragility, 
depending on the polymer-NP interaction, and that this 
effect becomes more pronounced with increasing (j). We 
relate these changes to fluctuations of {v?), which can be 
interpreted as a change in the fluctuations of the local 
molecular stiffness [H]. We also flnd that the Buchenau 
relation [15, holds with remarkable generality for all 
and interaction types, conflrming the relationship be- 
tween high frequency relaxation and structural relaxation 
under rather general circumstances. 

Our findings are based on equilibrium molecular dy- 
namics simulations of a nanoparticle surrounded by a 
dense polymer melt, as well as simulations of a pure melt 
for comparison purposes. We utilize periodic bound- 
ary conditions so that our results correspond to a uni- 
form dispersion of NP. The polymers are modeled by a 
well-studied bead-spring model |16j . All monomer pairs 
interact via a Lennard Jones (LJ) potential Vlj, and 
bonded monomers along a chain are connected via a 
FENE anharmonic spring potential. The NP consists of 
356 Lennard- Jones particles bonded to form an icosahe- 
dral NP; the facet size of the NP roughly equals the equi- 
librium end-to-end distance for a chain of 20 monomers. 
Details of the simulation protocol and our model poten- 
tials can be found in ref. [9]; further studies of the clus- 
tering and mechanical properties related to this model 
are presented in [TTJ [TS] . 

We simulate systems with 100, 200, or 400 chains of 
M = 20 monomers each (for totals of = 2000, 4000, 
and 8000 monomers) to address the effect of varying the 
NP volume fraction. Under constant pressure conditions. 
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FIG. 1: Relaxation time r as a function of T for each (f> nor- 
malized by Tpurc for the pure melt. The inset shows the raw 
data for r. Attractive interactions lead to increases in r, while 
non-attractive interactions decrease r. In both cases, the ef- 
fect is more pronounced with increasing concentration. The 
concentrations are (f) = 0.0426 (red Q s-nd orange A), 0.0817 
(green □ and brown <) and 0.151 (blue O and violet v)- The 
pure melt is indicated in black. The fits of the VFT form to 
the data deviate by at most 0.5 %. 

the addition of nanoparticles can give rise to a change 
in the overall melt density. A slight change in density 
can cause a significant but trivial change in the dynamic 
properties relative to the pure melt. In order to probe 
only changes caused by the interactions between the NP 
and the polymer melt, we have matched the density of 
monomers far from the NP with that of the pure polymer 
melt. 

To quantify changes in Tg and fragility, we evaluate ef- 
fect of (j) and the polymer-NP interactions on r, measured 
from the decay of the coherent intermediate scattering 
function (see supplementary information for a technical 
description). The effects of these interactions on r for 
some (j) were already presented in ref . [9] ; here we provide 
additional simulation data and expand the analysis con- 
sidered before to include the effect of the NP on fragility. 
As expected, Fig.[l]shows that attractive polymer-NP in- 
teractions slow the relaxation (r becomes larger), while 
non-attractive polymer-NP interactions give rise to an 
increased rate of relaxation (r becomes smaller). The ef- 
fect of (j) is more clearly seen by rescaling r by Tpure of 
the pure melt, which shows that r can be altered by a 
factor of more than an order of magnitude on cooling. 
The effect of the NP is more pronounced at low T, since 
at high T the polymer-NP interaction potential strength 
is weak in comparison with the system kinetic energy. 

We next examine how these changes in r affect Tg 
and fragility. The inset of Fig. |2] confirms that Tg in- 
creases when there is attraction, and decreases with non- 
attractive interactions. To estimate Tg, we fit the data 
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FIG. 2: Fragility dependence on </!> relative to the pure melt. 
We consider four different measures of fragility, which are dis- 
cussed in the text. All measures show the same qualitative 
trend: namely, the system with attractive polymer-NP in- 
teractions becomes more fragile, while the system with non- 
attractive interactions becomes less fragile (stronger). The 
inset shows the changes in Tg for the different (j) and interac- 
tions. 

using the Vogel-Fulcher-Tammann (VFT) expression |llj 

r = roe^/™-i). (1) 

To is an extrapolated divergence temperature of t, while 
D provides a measure of the fragility. We use the VFT 
fit to estimate Tg based on the condition that T{Tg) — 
100 s (the canonical definition of the laboratory glass 
transition [11 ), assuming the one time unit in standard 
LJ reduced units corresponds to I ps (reduced units are 
defined the in the supplementary information). 

Since Tg changes have been the focus of previous work, 
we emphasize how polymer-NP interactions and affect 
fragility. Since there is no single agreed upon measure of 
fragility, we consider several different measures. First, as 
indicated above, the parameter D from a VFT fit to r is 
widely utilized; specifically, a larger value of D indicates 
a stronger (less fragile) glass-forming fluid so that 1/D 
increases with increasing fragility. Another common esti- 
mate of the fragility is defined by the temperature ratios 
To/Tg or Tg/Tc- We estimate Tc using the power-law 
form T ~ (T/Tc — in an appropriate temperature 
range (see supplementary information for fitting details) . 
Since strong systems should have relatively weak tem- 
perature dependence of r approaching Tg, larger values 
of Tg/Tg or Tg/Tc correspond to more fragile systems. 
Finally, we use the VFT fit to estimate the fragility from 
the T dependence of r near Tg using the most commonly 
advocated fragihty definition [19 , 

m^d{lnT)/diTg/T)\T^. (2) 

For strong systems, the rate of change of r with respect 
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FIG. 3: T dependence of l/{u^) for different cj> normalized by 
I/(w^) for tlie pure melt. Note the parallelism to fig-[l] 

to T is smaller than that of fragile systems; hence m is 
larger for more fragile glass-forming fluids. We must rely 
on an extrapolation of the VFT fit to determine m, so 
that caution should be exercised interpreting the precise 
values of our m estimates. 

We summarize the results for the various fragility met- 
rics in Fig. [2] where we find that attractive polymer-NP 
interactions lead to more fragile glass formation as a func- 
tion of 0; conversely, non-attractive polymer-NP interac- 
tions lead to stronger glass formation. These changes 
are non-trivial, since Tg is independent or anticorrelated 
with fragility in some systems j20j . Our findings for the 
changes of fragility are consistent with several experimen- 
tal studies. Bansal et al. [7] found that dispersions of 
NP having repulsive interactions caused Tg to decrease, 
accompanied by an appreciable broadening of the glass 
transition region. Though they did not interpret their re- 
sults in terms of fragility, the increased breadth is indica- 
tive of increased strength (decreased fragility) , as we ob- 
serve. For the case of fuUerenes dispersed in polystyrene, 
Sanz et al. [H] reported behavior expected for attractive 
polymer-NP interactions, namely an increase in Tg, ac- 
companied by an increased fragility. Other studies have 
indicated no detectable or only a small change in the 
fragility with NP additive for small cj) [12] . The apparent 
absence of fragility changes at small <j> is also consistent 
with our results, which likewise show nearly undetectable 
changes in this concentration range. 

Having established the variation in Tg and fragility for 
the various systems, we now examine how these changes 
can be understood from the high frequency melt dynam- 
ics. Both experiments and simulations have shown that 
the Debye- Waller factor (m^) can be related to the low fre- 
quency relaxation associated with large scale structural 
relaxation |231l24j . (u^) measures monomer displacement 
on a time scale over which the particles are caged by their 
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FIG. 4: T dependence of the inverse variance l/cr associated 
with (u^) for different concentrations normalized by a for the 
pure melt. Symbols are the same as figs. [l] and [3] The changes 
in a mirror the changes of fragility, as argued in the text. The 
inset of the figure provides a reference for the T dependence of 
the distribution P{y?) for the system with attractive polymer- 
NP interactions. 

neighbors, and is accessible from both x-ray and neutron 
scattering measurements [25j . Since {v?) is usually deter- 
mined experimentally at a fixed instrumental time cor- 
responding to the time scale on the order of vibrational 
motion of the molecules, we determine the mean-squared 
chain segment displacement (MSD) at a vibrational time 
scale, specifically at a reduced time equal to 1.53, ~ the 
mean collision time. 

Fig. [3] shows that, when normalized by the behavior 
of the pure melt, the changes in (u^)purc/ (u^) are very 
similar to those observed for log(T/rpurc)- Evidently, the 
(u^) is sensitive not only to changes in 0, but also to 
the polymer-NP interactions. This affirms that (m^) is a 
potentially useful indicator of changes in low frequency 
relaxation. 

To understand how the changes in fragility relate to the 
high frequency dynamics, we must look beyond changes 
in (u^). It has been argued [T^ [55] that the changes 
in fragility can be connected with the efficiency of local 
packing, which is accompanied by changes in the fluctu- 
ations in the local moduli. Since the inverse of (u^) is a 
measure of the local liquid rigidity at high frequency [14] , 
we expect that fragility changes should be reflected in the 
variance cr^ = {u"^) — (u^)^ of the Debye- Waller factor. 

Figure [4] shows that, when normalized by the value 
for the pure melt, a is larger for stronger glass forma- 
tion (the non-attractive interactions), and smaller in the 
more fragile case (attractive interactions) in the T range 
considered. If we consider the normalized inverse of this 
quantity crpuj-o/c, we can examine to the relative variance 
in the local stiffness. These results show that the stronger 
systems have relatively smaller fluctuations in the local 
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FIG. 5; Parametric relation between r and {v?'). The inset 
shows the raw data, and the main plot shows the data collapse 
when the axes are scaled using the constants of the Buchenau 
relation. 

stiffness, while the more fragile systems have relatively 
larger fluctuations. This is consistent with the idea that 
the the stronger systems should be better packed, and 
vice- versa for the more fragile systems. Hence, l/cr ap- 
pears to be reflective of a change in fragility. 

Finally, we explore the proposed quantitative relation 
between r and (u^) provided by the "Buchenau relation," 



(3) 



where tb and Mq are system dependent constants. This 
relation has been verified for a number of systems, in- 
cluding the pure polymer melt currently under investiga- 
tion [21] . Larini et al. [23] have shown that a generaliza- 
tion of this relation seems to hold for diverse collection of 
simulated and real fluids over a wide range of tempera- 
tures, suggesting an amazing generality of the Buchenau 
expression, supporting the proposal that r is a universal 
function of {v?). The inset in Fig. [s] shows a striking 
agreement of our relaxation data with the Buchenau re- 
lation, the correlation becoming better at lower T. Devi- 
ations at high T can be expected since {v?) becomes pro- 
gressively ill-defined as T increases. More significantly. 



if we reduce r and 



by the corresponding fit param- 



eters Tb and Ug, we find that the data for all (j) and T 
can all be collapsed to a single master curve (main part 
of figure). This reinforces the somewhat surprising fact 
that the slow dynamics of the system are intimately con- 
nected with the high frequency vibrational properties of 
the system, as epitomized by {u^)- 

Similar to our findings, Riggleman et al. |26j have 
shown that changes in fragility associated with the addi- 
tion of molecular additives are linked to changes in the 
shear modulus of the material, another high frequency 
dynamical property. Additionally, Papakonstantopoulos 



et al. [57] have suggested a link between fragility varia- 
tions and elastic constant (shear modulus) fiuctuations in 
polymer nanocomposites in the glass state, but this was 
not verified through a direct determination of fragility. 

Care must be taken in comparing the glass below Tg 
with a liquid approaching Tg from above, as recent results 
suggest opposite trends may result. For example, Riggle- 
man et al. [26] show that a decrease of the fragility actu- 
ally results in a stiffening below Tg , but a softening above 
Tg, consistent with our results above Tg. Such an in- 
version between relative material stiffness and relaxation 
times near Tg has been suggested to be a general phe- 
nomenon [35], having potential applications in relation 
to scratch resistance of films, stabilizing nano-fabricated 
structures formed by lithography and nano-imprint tech- 
nology, and the preservation of drugs in glass-forming 
preservative formulations |28j . 
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